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Abstract 

An analytical theory for the efficiency of scattering-induced 
transitions from a random to a channeled state (feed-in) in bent 
crystals is derived. The predictions from the theory are in good 
agreement with experiment and Monte Carlo simulations. 

1 Introduction 

Steering of particle beams by means of channeling in bent crystals |1J is 
rapidly evolving technique for particle accelerators 0. The particle motion 
in a crystal is influenced by a series of collisions with the crystal constituents, 
which may cause particle transitions between channeled and random states. 
For any trajectory of a particle in a crystal a time-reversed trajectory is 
possible. The starting point of a trajectory becomes the final one, and vice 
versa. This leads to the idea of reversibility of transition processes [Q. In the 
depth of a crystal, besides the particles leaving the channeling mode (dechan- 
neling or feeding out) , there may be particles entering the channeling mode 
(respectively, feeding in, also known as volume capture). The mechanisms, 
responsible for these two opposite processes, are essentially the same. Here 
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we only discuss the scattering-induced transitions. Another mechanism of 
feeding- in, namely the centrifugal effects, is discussed elsewhere [Q. 

The feed-out process is well described by both the analytical theory and 
the computer simulations (Ref. |2j and references therein). The feed-in has 
been studied numerically, while the analytical theory was missing. Here we 
derive an explicit formula for the efficiency of the scattering-induced feeding- 
in in bent crystals. 



2 Feed-in efficiency in bent crystals 

Let us consider a beam with a uniform angular distribution, 1/2$, incident 
on a bent crystal with curvature 1/R. The fraction of beam channeled is a 
function of the crystal depth z: 

f(z) = 2 ^^A B ( P v/R)F D (z) (i) 

where As(pv/R) describes the reduction of the bent-crystal acceptance for 
a particle with momentum p and velocity v ; 8 C and x c are the critical angle 
and transverse position, respectively; d p is the interplanar spacing. The 
factor 7r/4 is exact in a harmonic case; for the realistic potential it should 
be replaced with «0.8. The Fd(z) describes the feed-out of the channeled 
particles due to scattering; usually one writes F D (z) as exp(— z/L D ) (L D is 
the dechanneling length), and so F' D (z) = —F D (z)/L D . 

The number of particles dechanneled over the length 5z equals 

df(z) 2x c tt6 c . . f . . 2x c tt6 c F d (z) , 

- — = ~ir;^ A *(pv/R)F D (z) = —--A B {po/R)— (2) 

The particles dechanneled over dz are exiting in the angular range dO = dz/R. 
Therefore the angular distribution downstream of the crystal is 

d f v d f( z ) 2x c tt6 c R 

d§ = = -d;^ A B(pv/R)j^F D (z) (3) 

We don't take into account a small extra spreading, ±# c , of the dechanneled 
particles. Therefore, for big R (when Lp/R is comparable to 9 C ), Eq. (|3|) 
overestimates the phase density. 
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Let us now consider the same beam incident on the same crystal in the 
reverse direction. Now the particles with the upstream parameters 
equal to the downstream parameters (xf,8f) of the dechanneled particles 
from the preceding captured along the same (reversed) trajectories. 

By consideration, the number of particles which have experienced transitions 
from the channeled states in the former case, is equal to the number of 
transitions to the channeled states in the latter case (as the trajectories are 
the same). 

Therefore the number of particles captured from the interval dO and then 
transmitted in the channeled states (over the length z) to the crystal face, 
is given by Eq. (|3|). We write this number as wsF D (z), the product of the 
capture probability and the transmission factor. The transmission factor 
Fjj(z), for particles channeled in the same states, is the same irrespective 
of the direction of motion. Normalizing Eq. (|3[) to the number of particles 
incident on the crystal in this angular range, d8/2Q, one obtains the capture 
probability: 

, df(z) 1 nx c R9 C t , , 
w s = 2$^pTryr = ~ , c /m A B (pv R) 4 
d9 F D (z) d p L D (pv/R) 

For a harmonic potential, Ab={1 — R c /R) 2 where R c is the critical radius; at 
the same time, Lb is reduced in a bent crystal by the same factor (1 — R C /R) 2 
relative to Lb in unbent crystal || . For a realistic potential the ratio of the 
two factors is ~1. Hence we can omit A B in (|j) and imply Lb value for a 
straight crystal; another simplification is x c /d p ml/2. One obtains: 

ttx c R9 C ^ ii R9 C ^ 

dp Lb 2 Lb 

The fact that ws should be of the order of R9 c /Lb has been earlier found 

from a simple qualitative consideration |J. To see the explicit dependence 

of Eq. on the properties of crystal and the particle energy, one can use 
the formulas for 9 C [|J 

\ pv J 
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and for Lb @ 

j _ 256 pv a TF d p 

D ~ 9vr 2 ' ln(2m e c 2 7 //) - 1 ' ^e 2 ' 1 j 

where N is the volume density of atoms, Z the atomic number, Z^e the 
particle charge, cltf the Thomas-Fermi screening distance, m e the electron 
rest mass, / ~ 16Z a9 eV is the ionization potential, 7 the particle Lorentz 
factor. Eq. (§) takes the form: 



^ = W J .fl.-^-(ln(2mecf» 7 /i)-l). (8) 



It should be mentioned that the length Lb for dechanneling from the 'stable 
states' is well defined and has been measured in many experiments, in good 
agreement with Eq. ([/]) ||. Therefore, the rate of feeding-in to the same 
states is an equally- well defined quantity. 

The feed-in rate has been studied in bent crystals experimentally at up 
to 70 GeV j7], It was found that ws is proportional to R and to p~ 3 / 2 , 
which is in perfect agreement with Eq. The Table 1 shows the prob- 

abilities of capture into the 'stable states' (which dechannel by exponential 
law) measured in the experiment at 70 GeV ||, found in the Monte Carlo 
simulation M, and calculated by Eq. (^). There is agreement within ~10-20 
%, i.e. within the errors of the experiment and simulation. 

Few remarks on Eq. (|5]). The feed-out rate at a very small depth in 

crystal is larger due to the 'short-lived' states with high transverse energies. 
Then we find from our consideration that the feed-in rate ~ R8 c /Ld to 
such states is respectively larger, which is qualitatively obvious. 

With R increase, Eq. ([|) may become > 1. The reason is the neglect of 
the ±9 C spreading in Eq. (§). If one sums quadratically the divergences, ±9 C 
and that from (|J), then the capture probability does not exceed 1. However, 
notice that ws is as small as ~1 % in the experimental practice. 

The reversibility relation (|5]) is valid irrespective of the mechanism of 
feeding-out (-in). In the presence of dislocations, Lb is small and there 
should be a high rate of feeding-in, in accordance with Eq. (0). This relation 



was observed indeed in the experiment |10| with dislocation-contaminated 



crystal of germanium and the 70 GeV proton beam. 
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Table 1: The probability (in %) of the feed-in into "stable states", for 70 GeV 
proton in Si crystal bent with R=3 m, from the experiment, simulation, and 
Eq. (5). 



Crystal 


Eq. (5) 


Simulation 


Experiment 


111 


0.20 


0.17±0.02 


0.23 


110 


0.22 


0.23±0.02 





3 Applications 

Although the feed-in processes cannot increase the beam bending efficiency, 
they contribute to a wide angular acceptance (^> 6 C ) crystal deflector. The 
bending efficiency of such deflector can be designed, and be varied in a broad 
range of values. Moreover, this efficiency is independent of the incident beam 
divergence (in contrast to the regular case of the entry-face capture). This 
may be valuable, when the beam attenuation, stability, and low background 
are important issues. Obviously, the secondary particles produced in colli- 
sions with crystal nuclei can be captured into the channeling mode (from 
the crystal bulk) through the feed-in processes only. This may be applicable, 
e.g., in some ideas of crystal application for experiments in particle physics 

El- 
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